Nonlinear dynamics in one dimension: 
On a criterion for coarsening and its temporal law 
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We develop a general criterion about coarsening for a class of nonlinear evolution equations de- 
scribing one dimensional pattern-forming systems. This criterion allows one to discriminate between 
the situation where a coarsening process takes place and the one where the wavelength is fixed in the 
course of time. An intermediate scenario may occur, namely 'interrupted coarsening'. The power 
of the criterion on which a brief account has been given [P. Politi and C. Misbah, Phys. Rev. Lett. 
92 , 090601 (2004)], and which we extend here to more general equations, lies in the fact that the 
statement about the occurrence of coarsening, or selection of a length scale, can be made by only 
inspecting the behavior of the branch of steady state periodic solutions. The criterion states that 
coarsening occurs if X'(A) > while a length scale selection prevails if A' (A) < 0, where A is the 
wavelength of the pattern and A is the amplitude of the profile (prime refers to differentiation) . This 
criterion is established thanks to the analysis of the phase diffusion equation of the pattern. We 
connect the phase diffusion coefficient D(X) (which carries a kinetic information) to A' (A), which 
refers to a pure steady state property. The relationship between kinetics and the behavior of the 
branch of steady state solutions is established fully analytically for several classes of equations. An- 
other important and new result which emerges here is that the exploitation of the phase diffusion 
coefficient enables us to determine in a rather straightforward manner the dynamical coarsening 
exponent. Our calculation, based on the idea that |D(A)| ~ A 2 /t, is exemplified on several nonlinear 
equations, showing that the exact exponent is captured. We are not aware of another method that 
so systematically provides the coarsening exponent. Contrary to many situations where the one 
dimensional character has proven essential for the derivation of the coarsening exponent, this idea 
can be used, in principle, at any dimension. Some speculations about the extension of the present 
results are outlined. 

PACS numbers: 05.70.Ln, 05.45.-a, 82.40.Ck, 02.30.Jr 



I. INTRODUCTION 

Pattern formation is ubiquitous in nature, and espe- 
cially for systems which are brought away from equi- 
librium. Examples are encountered in hydrodynam- 
ics, reaction-diffusion systems, interfacial problems, and 
so on. There is now an abundant literature on this 
topic 0, Q]. Generically, the first stage of pattern for- 
mation is the loss of stability of the homogeneous solu- 
tion against a spatially periodic modulation. This gen- 
erally occurs at a critical value of a control parameter, 
fi = /i c (where /i stands for the control parameter) and 
at a critical wavenumber q — q c . The dispersion rela- 
tion about the homogeneous solution (where perturba- 
tions are sought as e i i x + UJt ^ j n the vicinity of the critical 
point assumes, in most of pattern-forming systems, the 
following parabolic form (Fig. inset) 

LJ = S - (q - q c f (1) 
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where S is proportional to (/x — /i c ). For S < 0, u> < for 
all q's and the homogeneous state is stable. Conversely, 
for 5 > there is a band of wavevectors Aq = (q — q c ) = 
±V5 corresponding to unstable modes (Fig. 1), so that 
infinitesimal perturbations grow exponentially with time 
until nonlinear effects can no longer be ignored. In the 
vicinity of the bifurcation point (6 = 0) only the princi- 
pal harmonic with q = q c is unstable, while all other har- 
monics are stable. For example, Rayleigh-Benard convec- 
tion, Turing systems, and so on, fall within this category, 
and their nonlinear evolution equation is universal in the 
vicinity of the bifurcation point. If the field of interest 
(say a chemical concentration) is written as A(x,t)e lqcX , 
where A is a complex slowly varying amplitude, then A 
obeys the canonical equation 

d t A = A + d xx A- \A\ 2 A (2) 

where it is supposed that the coefficient of the cubic term 
is negative to ensure a nonlinear saturation. Because the 
band of active modes is narrow and centered around the 
principal harmonic, no coarsening can occur, and the pat- 
tern will select a given length, which is often close to that 
of the linearly fastest growing mode. However, the ampli- 
tude equation above exhibits a phase instability, known 
under the Eckhaus instability Q] , stating that among the 
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band of allowed states, |Ag| = s/5, only those modes 
whose wavevectors satisfy \Aq\ < arc stable with 

respect to a wavelength modulation. 

There arc many other situations where the bifurca- 
tion wavenumber q c —* and therefore a separation of 
a slow amplitude and a fast oscillation is illegitimate. 
Contray to the case Q, where the field can be written 
as A(x : t)e lqaX with A being supposed to vary slowly in 
space and time, if q c — > the supposed fast oscillation, 
e lqaX , becomes slow as well and a separation of A does 
not make a sense anymore. In this case, a generic form 
of the dispersion relation is (Fig. main) 

uj = dq 2 - q 4 . (3) 

A third situation is the one where the dispersion relation 
takes the form (Fig. [21 inset) 

Lu = 5-q 2 . (4) 

In both cases, Eqs. (l.'il 111 , the instability occurs for S > 0, 
and the band of unstable modes extends from q = to 
q = V5. That is to mean, there is an infinite number of 
unstable harmonics: if q is the wavenumber of an unsta- 
ble mode, then also q/2,q/3, . . . are unstable. Examples 
that fall in this category are numerous 3]: the air-liquid 
interface in a thin film falling on an inclined plane, flame 
fronts, step dynamics in step flow growth, sand ripples, 
and so on, or simply the Ginzburg-Landau equation J5J, 
which corresponds to case (@J). 

Dispersion relation has an extra factor of q 2 which 
is often due to a conservation constraint (see also below). 
Because of the dispersion form, constant plus a quadratic 
term, Q might formally resemble (Q. However, an im- 
portant caution must be taken: in Q it must be remem- 
bered that q should remain close to q Cl so that only one 
harmonic is active, while in n o such a restriction is 
made, and therefore q can be as close as possible to zero, 
leading to a highly nonlinear dynamics. 

Other types of dispersion relations which may arise, 
and which are worth of mention, correspond to situa- 
tions where w = 5\q\ — q 2 , or id = S\q\ — \q\ 3 , leading 
also to a vigorous mode mixing, for the same reasons 
evoked above. The occurrence of a non analytic dis- 
persion relation with \q\ is a consequence of long range 
interactions If the unstable band extends down to 
q — 0, the appropriate form of the evolution equation is 
not an amplitude equation for a slowly varying function 
A, but rather a partial nonlinear differential equation, or 
an integro-differential equation, for the full field of inter- 
est, say h(x,t) if one has in mind a front profile. 

A prominent example of a partial differential equation 
is the Kuramoto-Sivashinsky |g,|g,|2j (KS) equation 

dtu = —d xx u — d xxxx u + ud x u (5) 



In case J1J an equation similar to J2J may arise, but it describes 
the full field and not just the envelope. 
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Swift-Hohenberg eq. 
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FIG. 1: Inset: Dispersion curve JTJ: u = 8- (q-q c ) 2 . There 
is loss of stability at a critical wavenumber q c , when 5 = 
(dotted line). For 8 > 0, the unstable band extends from 
q — qi to q — q2 (full line). With increasing 5 the unstable 
region widens (dashed line). The parabolic shape of w(q) is 
an approximation, valid close to its maximum. This applies, 
e.g., to the dispersion curve 1381 of the Swift-Hohenberg eq. 
(see the main figure): w = 5 — (1 — q 2 ) 2 [q e = 1]. When 8 — 1 
(dashed line) the unstable band extends down to q = and 
uj(q) resembles the dispersion curve of the Cahn-Hilliard eq. 
sec © and Fig.0. 
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FIG. 2: Main: dispersion curve ©, valid, e.g., for the Cahn- 
Hilliard equation: u — 8q 2 — q 4 . Dotted line: below threshold 
(8 < 0). Full line: just above threshold. Dashed line: well 
above threshold. The vanishing of w(0) for any 8 is a con- 
sequence of the traslational invariance of the CH eq. in the 
'growth' direction. Inset: The Ginzburg-Landau eq. is one 
case where such invariance is absent and the dispersion curve 
has the form Q: to = 8 — q 2 . 



which leads to spatio-themporal chaos. Note that by set- 
ting u — d x h we obtain an equivalent form of this equa- 
tion, namely d t h = -d xx h - d xxxx h + (d x h) 2 /2. This 
equation arises in several contexts: liquid films flowing 
down an inclined plane ||, flame fronts 0, step flow 
growth of a crystal surface |3] ■ 

Complex dynamics such as chaos, coarsening, etc . . . , 
are naturally expected if modes of arbitrarily large wave- 
length are unstable. However, these dynamics may occur 
for systems characterized by the dispersion relation Q 
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as well, if the system is further driven away from the 
critical point (i.e., if q 2 ^S> qi, see Fig.QJ because higher 
and higher harmonics become active. We may expect, 
for example, coarsening to become possible up to a total 
wavelength of the order of 2iz/q\. 

For systems which are at global equilibrium the nonlin- 
earity ud x u is not allowed, and a prototypical equation 
having the dispersion relation © is the Cahn-Hilliard 
equation 

d t u = -d xx [u + d xx u - u 3 ]. (6) 

The linear terms are identical to the KS one, and the dif- 
ference arises from the nonlinear term. Note that if dy- 
namics is not subject to a conservation constraint, (— d xx ) 
on the right hand side is absent, and the dispersion rela- 
tion is given by Eq. |0J. The resulting equation is given 
by (01 for a real A and it is called real Ginzburg-Landau 
(GL) equation or Allcn-Cahn equation. 

The KS equation, or its conserved form (obtained by 
applying d xx on the right hand side), was suspected for 
a long time to arise as the generic nonlinear evolution 
equation for noncquilibrium systems (the quadratic term 
is non variational in that it can not be written as a func- 
tional derivative) whenever a dispersion relation is of type 
©. Several recent studies, especially in Molecular Beam 
Epitaxy (MBE) , have revealed an increasing evidence for 
the occurrence of completely new types of equations, with 
a variety of dynamics: besides chaos, there are ordered 
multisoliton [& || solutions, coarsening ^(j> freezing of 
the wavelength accompanied by a perpetual increase of 
the amplitude [llj. Moreover, eq uations bearing strong 
resemblance with each other |l2( exhibit a completely 
different dynamics. Thus it is highly desirable to extract 
some general criteria that allow one to discriminate be- 
tween various dynamics. 

A central question that has remained open so far, 
and which has been the subject of a recent brief expo- 
sition was the understanding of the general condi- 
tions under which dynamics should lead to coarsening, 
or rather to a selection of a length scale. In this paper 
we shall generalize our proof presented in [T^ | to a larger 
number of classes of nonlinear equations, for which the 
same general criterion applies: the sign of the phase dif- 
fusion coefficient D is linked to a property of the steady 
state branch. More precisely, the sign of D is shown to 
be the opposite of the sign of X'(A), the derivative of the 
wavelength A of the steady state with respect its ampli- 
tude A. Therefore, coarsening occurs if (and only if) the 
wavelength increases with the amplitude. 

Another important new feature that constitutes a sub- 
ject of this paper, is the fact that the exploitation of the 
phase diffusion coefficient D(X) will allow us to derive an- 
alytically the coarsening exponent, i.e. the law according 
to which the wavelength of the pattern increases in time. 
For all known nonlinear equations whose dispersion rela- 
tion has the form © or 10} and display coarsening, we 
have obtained the exact value of the coarsening expo- 
nent, and we predict exponents for other non exploited 



yet equations. An important point is that this is expected 
to work at any dimension. Indeed, the derivation of the 
phase equation can be done in higher dimension as well. 
If our criterion, based on the idea that |-D(A)| ^ \ 2 /t, 
remains valid at higher dimensions, it should become a 
precious tool for a straightforward derivation of the coars- 
ening exponent at any dimension. 



II. THE PHASE EQUATION METHOD 

A. Generality 

Coarsening of an ordered pattern occurs if steady state 
periodic solutions are unstable with respect to wave- 
length fluctuations. The phase equation method al- 
lows to study in a perturbative way the modulations of 
the phase <f> of the pattern. For a periodic structure of 
period A, (f> = qx, where q = 2ir/\ is a constant. If wc 
perturb this structure, q acquires a space and time de- 
pendence and the phase <p is seen to satisfy a diffusion 
equation, dt<j) = Dd xx (f>. The quantity D, called phase 
diffusion coefficient, is a function of the steady state so- 
lutions and its sign determines the stable (D > 0) or 
unstable (D < 0) character of a wavelength perturba- 
tion. 

A negative value of D induces a coarsening process, 2 
whose typical time and length scales are related by 
|-D(A)| ~ A 2 /i, as simply derived from the solution of 
the phase diffusion equation: this relation allows to find 
the coarsening law X(t). Therefore, the phase equation 
method not only allows to determine if certain classes 
of partial differential equations (PDE) display coarsen- 
ing or not; it also allows to find the coarsening laws, 
when D < 0. In the rest of this section, we are going 
to offer a short exposition of the phase equation method 
without referring to any specific PDE. Explicit evolution 
equations will be treated in the next sections, with some 
calculations relegated to the appendix. 

Let us consider a general PDE of the form 3 

d t u(x,t) =J\f[u] (7) 

where Af is an unspecified nonlinear operator, which is 
assumed not to depend explicitly on space and time. 
uo(x) is a periodic steady state solution: A/"[mo] = and 
u (x + A) = u (x). 

When studying the perturbation of a steady state, it is 
useful to separate a fast spatial variable from slow time 
and space dependencies. The stationary solution uq does 



2 In principle a negative D could entail also a decreasing of A 
(splitting). However, this is inconsistent with the result dT/d\ < 
(see Sec.|VJ and with the stability of the fiat interface at small 
length scales. 

3 Coarsening scenarios are not affected by the presence of noise, 
which is not taken into account throughout the article. 
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not depend on time and it has a fast spatial dependence, 
which is conveniently expressed through the phase <fi = 
qx. Once we perturb the stationary solution, 

u = u + eu% + . . . , (8) 

the wavevector q = d x (j> gets a slow space and time de- 
pendence: q = q(X,T), where X = ex and T = e a t. 
Because of the diffusive character of the phase variable, 
the exponent a is equal to two. Space and time deriva- 
tives now read 

d x = qd^ + ed x (9a) 
d t = e(d T ^ (9b) 

where the second order term in the latter equation (e 2 dr) 
has been neglected. Finally, along with the phase <j> it is 
useful to introduce the slow phase i/j(X, T) = ecf>(x, t), so 
that q = dxip- 

Replacing the u— expansion JSJ and the derivates © 
with respect to the new variables in Eq. J7J), we find an 
e— expansion which must be vanished term by term. The 
zero order equation is trivial, A/o[uo] = 0: this equation is 
just the rephrasing of the time-independent equation in 
terms of the phase variable <fi (the subscript in Afo means 
that Eqs. 10 have been applied at zero order in e, i.e. 
dx = qdj,). 

The first order equation is more complicated, because 
both the operator AT and the solution u are e— expanded. 
On very general grounds, we can rewrite dtu(x, t) = Af[u] 
as 

e(d T i))d 4> u (s = (Af + eA/i)[u + em] (10) 

where Afi comes from first order contributions to the 
derivatives JjJJl- If we use the Frechet derivative 0, Co, 
defined through the relation 

Af [u + em] = Af [u ] + e£ [«i] + 0{e 2 ) (11) 

we get 

C~o[ux] = (dr^d^uo - Mi[u ] = g(u ,q,ijj). (12) 

At first order, therefore, we get an heterogeneous linear 
equation (the Frechet derivative of a nonlinear operator 
is linear). The translational invariance of the operator 
Af guarantees that d^u^ is solution of the homogeneous 
equation: according to the Fredholm alternative theo- 
rem [lGl . a solution for the heterogeneous equation may 
exist only if g is orthogonal to the null space of the ad- 
joint operator £' . In simple words, if C'[v] = 0, v and g 
must be orthogonal. This condition, see Eq. i|12|) . reads 

(v, d^uoidr^ = (v,A7i[uo\), (13) 

where 4 (f,g) - {2^ d^f*g. 



4 Sometimes we may also write {/(</>)) to mean (27r) 1 $ d(f>f '. 



It happens that A/i[uo] is proportional to dxq = dxxi 3 ) 
and the previous equation has the form of a diffusion 
equation for the phase tp, 

d t iP = Ddxx^- (14) 

B. Applications 

1. The generalized Ginzburg- Landau equation 

The (real) Ginzburg-Landau equation is written as 

d t u = (u - u 3 ) + u xx (15) 

whose linear spectrum, for an excitation u(x,t) = 
exp(wt + iqx) is uj(q) = 1 — q 2 . This equation is the 
prototype for the evolution of a nonconserved order pa- 
rameter with two equivalent stable solutions, u = ±1. 
Starting from the trivial solution u = 0, we have a linear 
instability leading to a logarithmically slow coarsening 
process |l7j . 

This equation can be easily generalized to 

d t u = B{u) + G{u) u xx , (16) 

which will therefore be called generalized Ginzburg- 
Landau (g'GL) equation. If B(u) ~ u and G(u) ~ 1 
for small u, the linear spectrum is unmodified, but 
the nonlinear behavior can be totally different, de- 
pending on the full (positive) expressions of B(u) and 
G(u). Steady states are determined by the relation 
u'q(x) = —B(u )/G(u ), so they correspond to the tra- 
jectories of a classical particle moving under the force 
-B(uo)/G(u ). 

Now, let us apply the expansions © for the order pa- 
rameter and © for the derivatives to Eq. l(TB|) . The first 
and second spatial derivatives can also be written as 

d x = qdcf, + eipxxd q (17a) 
d xx = q 2 d ct>ct> + e'ifjxx(2qd q + l)d < p (17b) 

As anticipated in the previous section, the zero and 
first order equations read A/q[uq] = and £o[ui] = .<?, 
where 

Af [uo] = B(u Q ) + G(u Q )q 2 d lf>< puo (18) 

is the nonlinear operator defining the gGL equation, 

£oM = [B'(u ) + G'{u )q 2 (d u u ) + G{u Q )q 2 d <i><j> }u 1 

(19) 

is its Frechet derivative, and 

g = {dr^uo - (dxx^)G{u ){2qd q + l)d^u . (20) 

Because of translational invariance, £o[^uo] = 0. Its 
adjoint is easily found to be 

4M = q 2 d^[vG(u )] + [B'(uq) + G'{u )q 2 (d^u )]v. 

(21) 
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If we define w = vG(uq), the equation Cq[v] = is iden- 
tical to £o[H = 0, so that we can choose w = d^uo and 
v = d^uo/Giuo). 

The orthogonality condition between v and g reads 

(dTi/})(v,d^uo) - (dxxip)(v,G(u )(2qd q + l)cfyuo) = 

(22) 

and replacing the explicit expression for v, we get the 
phase diffusion equation 



with 



D = 



\ G(«o) 



El 

Do' 



(23) 



(24) 



Assuming a positive G, the sign of D is fixed by the 
increasing or decreasing character of (qid^uo) 2 ) with the 
wavevector q. Reversing to the old variable x, 



(?(^ ) 2 ) = i [ dx(u' ) 2 = f 
2tt J q 2tt 



(25) 



where J is the well known action variable, whose deriva- 
tive with respect to the 'energy' of the particle gives the 
period A. The following relations are easily established: 



l_£U_ 
27 ~dq~ 



X 3 



4tt 2 \dE 



8X 



A 3 B(A) fax 



4tt 2 G{A) \dA 



(26) 

where A is the amplitude of the oscillation, i.e. the (pos- 
itive) maximal value for uq(x). 

If G(u) = 1, a compact formula for D is 



D = - 



X 2 B{A) 
J(d A X) 



G{A) = I. 



(27) 



In conclusion, a coarsening process occurs (D < 0) 
if the wavelength of the steady states increases with in- 
creasing their amplitude. In App.^we make some gen- 
eral remarks on the behavior of A (^4) for several different 
potentials. 



2. The generalized Cahn-Hilliard equation 

The Cahn-Hilliard equation is the conserved 'version' 
of the Ginzburg-Landau equation, 



The linear regime corresponds to an exponential unsta- 
ble growth of such mode, with a rate u(q u ), followed by 
a logarithmic coarsening. 5 

The above equation can be made of wider application 
by considering the following generalized Cahn-Hilliard 
(gCH) equation 



d t u = -C(u)d xx [B(u) + G{u)u xx ]. 



(29) 



In Sec. lIIIDl we will discuss thoroughly the coarsening 
of this class of models, because of its relevance for the 
crystal growth of vicinal surfaces. 6 In that case, the local 
height z(x,t) of the steps satisfies the equation 



8 t z = -d x {B(m) + G(m)d x [C(m)d x m}} 



(30) 



where m = d x z. If we pass to the new variable 
u(m) = J Q dsC(s) and take the spatial derivative of 
the above equation, we get the gCH equation (|2*9"jl . It 
is worthnoting that steady states are given by the equa- 
tion B(uq) + G(uq)u'q = jo i where j is a constant de- 
termined by the condition (uq) = tuq that imposes the 
(conserved) average value of the order parameter. If 
steps are oriented along a high-symmetry orientation, 
too = = jo. In the following we are considering this 
case only, so the equation determining steady states, 
B(uq) + G(uo)u'q = 0, is the same as for the gGL equa- 
tion. 

If we proceed along the lines explained in Sec III Al and 
keep in mind notations used in Sec lIIB H the first order 
equation in the small parameter e reads 



- q 2 d tj><j> C Q [ui} = (drip) 



(31) 



C(u ) 

+ (dxx^)q 2 d<t><t>[G(u )(2qd q + l)fyuo]. 

According to the Fredholm alternative theorem, the 
right hand side must be orthogonal to the solution v of 
the equation 



(9^£ ) t H = C-ld^v = 0. 



(32) 



According to the results of the previous section, we 
know that d^v — (d r f > u )/G(u )- The orthogonality con- 
dition now reads 



(33) 



q 2 {v,d^[G(u )(2qd q + l)fyu ])(ft«-VO = 0. 



d t u = -d xx [(u -it 3 ) +u xx \. (28) 

The spatial average of the order parameter is time 
independent, d(u)/dt = 0, and the linear spectrum is 
oj{q) = q 2 — q 4 : it therefore has a maximum at a finite 
value q u = l/v2) called the most unstable wavevector. 



5 The coarsening of the nonconserved (Ginzburg-Landau) and con- 
served (Cahn-Hilliard) models differ if noise is present: \(t) ~ 
t 1 / 2 in the former case and \(t) ~ t 1 / 3 in the latter case. 

6 A vicinal surface is a surface which is slightly miscut from a high- 
symmetry orientation. It looks like a flight of stairs with steps 
of atomic height separating large terraces. 
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The quantity multiplying (dxxip) can be rewritten as 

(v,d(M>[G(uo)(2qd q + l)a^uo) = 
'' " n -,G(u Q )(2qdg + l)d^u Q ) = 



d q (q(d^U Q ) 2 ) (34) 
so we finally have ($tV0 = D(dxxi } ) with 



D = - 



q 2 d q (q(d 4> uo) 2 ) _ q 2 D 1 



C(u ) 



Do 



(35) 



In App.[B]we prove that the denominator D 2 is always 
positive. If C and G are (positive) constants the proof is 
straightforward, because —(vO^uq) = ((d^v)^) — (uq). 
The diffusion coefficient l|35|) for the gCH equation is 
therefore similar to the diffusion coefficient (|24H for the 
nonconserved gGL equation: their sign is determined by 
the increasing or decreasing character of X(A), the wave- 
length of the steady state, with respect to its amplitude. 
The q 2 term in the numerator of (|35|l is evidence of the 
conservation law, i.e., of the second derivative d xx in 
Eq. (|29[l . The denominators D 2 and D 2 differ: this is 
irrelevant for the sign of D, but it is relevant for the 
coarsening law. 

If C(u) = G(u) = 1, formulas simplify: D\ = 
-\ 3 B(A)/4n 2 {d A \) and D 2 = (u 2 ) = where I = 
§ Uq(x) has the same role as J = §{d x uo) 2 in the non- 
conserved model. Putting everything together we obtain 



D 



\ 2 B{A) 



C{A) = G(A) = 1. 



(36) 



3. The generalized Swift-Hohenberg equation 

Both the (generalized) GL and CH equations have a 
linear spectrum whose unstable band extends to q = 0, so 
that active modes of arbitrarily large wavelength exist. In 
the Swift- Hohcnbcrg equation we can tune a parameter 8 
so that to change the unstable band. The standard form 
of the equation is 



dtu = Su — (1 + d xx ) U — u 

= -d*u - 2d 2 u - (1 - S)u - u 3 



(37) 



Linear stability analysis for a single harmonic u(x, t) 
exp(ujt + iqx) gives the spectrum 



2\2 



tu(q) = 6 - (I - q 2 ) 



(38) 



and for positive S there is a finite unstable band (u>(q) > 
0) extending from q\ — VT— Vd to q 2 = vl + \/~8. For 
5=1, the unstable band is the interval (0, v2)- The most 
unstable wavevector is q u = 1 for any 5. For small 5 the 
unstable band is narrow; in fact, for 5 < 0.36, q± > q 2 /2 



and period doubling is not allowed. In other words study- 
ing coarsening for the Swift-Hohenberg equation close to 
the threshold 6 = is not very interesting: nonetheless 
we will write the phase diffusion equation for any S and 
for a generalized form of the Swift-Hohenberg equation 
as well. 

The zero order equation is easy to write 

A/o [uo] = ~q%u - 2q 2 d 2 u - (1 - 5)u - u 3 = (39) 

and the first order equation has the expected form 
£o[ui] = g, where 

Co = ~q% - 2q 2 8 2 - (1 - S) - 3u% (40) 

is the Frechet derivative of No and 

g = {d T ^)d^uo (41) 
+2(dxK^mq fl d q + 3q 2 )dlu Q + {2qd q + l)fyu„]. 

The operator Co is self-adjoint, so the solution of the 
homogeneous equation Cq[v] — is immediately found, 
because of the translational invariance of No along x: 
v = d^uo- We therefore have 



(cvvoa^o) 2 ) = 

-{dxx^)[{d^uo{Aq 3 d q + 6q 2 )dlu ) 
+2(d^u (2qd q + l)d,pu )}. 



(42) 



It is easy to check that both terms appearing in square 
brackets on the right hand side can be written as d q (. . .): 



(cVo(V3 g + 6q 2 )d$u ) = -3,(2^(^7/0)') (43a) 
{d 4> uo{2qd q + l)fyuo>] = dMd^uof) (43b) 

so that the phase diffusion coefficients reads 

d q [2q 3 ({d 2 uo) 2 )-2q{{d^Uo) 2 )] 



D = 



(44) 



{{d^uo) 2 ) 

Now let us generalize this result to the equation 

d t u = c nd> + P(u) (45) 



n>0 



where n is assumed even and P(u) an odd function, in 
order to preserve the symmetries x — > —x and u — * — u. 
Therefore n = 2k, with k > 1. We report here only the 
final result for the phase diffusion coefficient: 



D = 



1 



2((cVo) 2 



Y.i-lf^ncnd^id^uof) (46) 



The standard Swift-Hohenberg equation corresponds 
to c 2 = -2 and c 4 = 1. The quantity (-l) k+x nc n q n - 1 /2 
therefore gives 2q 3 for n = 4 ans — 2q for n = 2, as shown 
by Eq. flg). 
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III. THE COARSENING EXPONENT 

We now want to use the results obtained in the previ- 
ous section for the phase diffusion coefficient D in order 
to get the coarsening law X(t). In one dimensional sys- 
tems, noise may be relevant and change the coarsening 
law. In the following we will restrict our analysis to the 
deterministic equations. 

A negative D implies an unstable behavior of the phase 
diffusion equation, d t ip = — \D\d xx tlj, which displays an 
exponential growth (we have reversed to the old coordi- 
nates for the sake of clarity): ip = exp(i/r) exp(27ria;/A), 
with (27r) 2 |Z?| = A 2 /r (in the following the time scale r 
will just be written as t). The relation |-D(A)| pa A 2 /i will 
therefore be used to obtain the coarsening law \(t): it 
will be done for several models displaying the scenario of 
perpetual coarsening (i.e., A' (A) > for diverging A). 



A. The standard Ginzburg-Landau and 
Cahn-Hilliard models 



In order to apply Eq. 1(27(1 and find D for the noncon- 
served model, we also need B(A) = A — A 3 ps 2exp(— A) 
and J; the latter quantity is the classical action and it is 
easy to understand that it keeps finite in the limit A — > 1, 
because it is the area in the phase space enclosed in the 
limiting trajectory. So, the relation |-D(A)| ps X 2 /t gives 
d A X/B(A) pa t and finally 

A(i)sslni standard GL eq. (47) 

If we compare Eq. ijTH) with Eq. we find that 

the phase diffusion coefficient D for the nonconserved 
model is equal to that for the conserved model, once 
J = §{d x u^) 2 has been replaced by / = §(uq) 2 , which 
is of order A because uq is almost everywhere equal to 
one. This replacement is irrelevant in our case, because 
the relation exp(2A) ~ t, giving rise to Eq. (|4*7)l in the 
nonconserved case, is replaced by Aexp(2A) ~ t in the 
present, conserved case, so that a logarithmic coarsening 
is found again: 

X(t) pa lni standard CH eq. (48) 



It is well known [TtJ that in the absence of noise, both 
the nonconserved GL equation 1(151) and the conserved 
CH equation display logarithmic coarsening, X(t) pa 
Int. Let us remind that steady states correspond to the 
trajectories of a classical particle moving in the potential 
V{u) — u 2 /2 — u 4 /4. The wavelength of the steady state, 
i.e. the oscillation period, diverges as the amplitude A 
goes to one. This limit corresponds to the 'late stage' 
regime in the dynamical problem, and the profile of the 
order parameter is a sequence of kinks and antikinks. The 
kink (antikink) is the stationary solution u + (x) (u_(x)) 
which connects u = —1 (u = 1) at x = — oo to u = 1 (u = 
— 1) at x = oo, u±(x) = ± tanh(x/\/2). As coarsening 
proceeds, kinks and antikinks annihilate. 

It is convenient to introduce the new variable Q(x) — 
1 — u(x). In the mechanical analogy, if the particle leaves 
(with zero velocity) from u(0) = A = I — Qq, the quantity 
Q(x) = l — u(x) starts to grow exponentially (this follows 
from the tanh(a;/\/2) solution, 7 Q(x) ~ Qo exp(\/2a;)). 
The particle passes a diverging time close to x = ±1 be- 
cause x = ±1 correspond to the potential maxima where 
the velocity of the particle vanishes while its position re- 
main finite. Consequently, it is straightforward to write 8 
that A ~ ln(l/Qo) ~ — ln(l — A) and calculate the deriva- 
tive 8aX pa (1 — A) -1 pa exp A. 



7 This may also be directly seen by expanding the differential equa- 
tion d xx u + u — u s = about u = 1. Expansion to leading order 
in Q = 1 — u yields d xx Q — 2Q = 0, from which the solution 
Q = Q e^ x follows. 

8 It is possible to calculate more rigorously the period of oscilla- 
tion from the exact relation A = 2 J_a du/(d x u), which can be 
expressed in terms of an elliptic integral. In the limit A — * 1 it 
diverges logarithmically as — ln(l — A). 



B. Other models with uniform stable solutions 

The only features of the standard GL and CH mod- 
els which determine the coarsening laws (|47|l and 148(1 
are the presence of symmetric maxima in V(u) for fi- 
nite u = u = 1 and the quadratic behavior close to 
the maxima, V(u) pa V{u) — a(u — u) 2 . For any gen- 
eralization of such equations, d t u — B(u) + u xx and 
dtu = —d xx [B(u) + u xx ], if y{ u ) — J duB{u) still has 
such properties, the coarsening will be logarithmic in- 
dependently of the details of B(u). Since a quadratic 
behavior close to a maximum is general, we can conclude 
that a logarithmic coarsening is common to most of non- 
conserved and conserved models where B{u) vanishes for 
finite u. 

In this section we investigate situations where the as- 
sumption that V(u) has a quadratic behavior close to 
the maxima is relaxed. They corresponds to very special 
situations, nonetheless they are interesting in principle 
because X(t) is no more logarithmic and because these 
models constitute a bridge between the standard GL/CH 
models and the models without uniform stable solutions, 
discussed in the next section. Let us now consider a class 
of models where V(u) has a minimum in u — (linearly 
unstable profile) and two symmetric maxima in u = ±1, 
with V(u) ~ V(l) - a\u - If ((3 > 2) close to u = 1. 
Again, we use Q{x) = 1 — u(x) and Qq = 1 — A, A being 
the amplitude of the oscillation. 

The wavelength is given by 



A 



dQ 



1 



\/Q' 3 - Qg Q6 



and it is easily found that J 
B(A)~Q$-\ 



1, / ~ A, d A X 



(49) 



-0/2 
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1. The nonconserved case 



The relation 



gives 



* ' ' J{d A \) 



0-1 



(50) 



(51) 



Q'o 



so that X(t) ~ t n with n=((3- 2)/(3/3 - 2). 



2. The conserved case 

If the order parameter is conserved, we simply need to 
replace J ~ 1 with / ~ A in Eq. (|50|) , so as to obtain 



3/3-2 

t ~ AAt^ 



4/3-4 

At^ 



(52) 



The coarsening exponent is therefore equal to n = 
{13 - 2)/ (4/3 - 4). We observe that in the limit j3 -> 2 we 
recover logarithmic coarsening (n = 0) both in the non- 
conserved and conserved case, as it should be. We also 
remark that in the opposite limit (3 — ► oo we get n = 1/3 
in the nonconserved case and n = 1/4 in the conserved 
case, which make a bridge towards the models discussed 
in the next section. 



C. Models without uniform stable solutions 

The models considered in the previous subsection have 
uniform stable solutions, u = ±1: the linear instabil- 
ity of the trivial solution u = leads to the formation 
of domains where the order parameter is alternatively 
equal to ±1, separated by domain walls, called kinks 
and antikinks. This property is related to the fact that 
B(u) — u — u 3 vanishes for finite u (up to a sign, B(u) is 
the force in the mechanical analogy for the steady states). 

In the following we are considering a modified class of 
models, where B(u) vanishes in the limit u — > oo only, 
so that the potential V(u) = J duB{u) does not have 
maxima at finite u. Therefore, it is not possible to define 
'domains' wherein the order parameter takes a constant 
value. These models [l^, which may be relevant for the 
epitaxial growth of a high-symmetry crystal surface 01 > 
are defined as follows (a > 1): 



d t u = B{u) + u xx 
d t u = -d xx [B(u) + u x 
B(u) = 



V(u) 



(l+u 2 ) a 

1 1 
~2(a- 1) (l + u 2 )"" 1 



(53) 
(54) 

(55) 
(56) 



Steady states correspond to periodic oscillations of a 
particle in the potential V(u) around the minimum in 
u = 0. With increasing the amplitude A, the energy 
of the particle goes to zero as E ~ — yj~( 2a ~ 2 ) and the 
motion can be split into two parts: the motion in a fi- 
nite region around u — 0, which does not depend on A, 
and the motion for large u, where V(u) ~ l/u^ 2a ^ and 
simple dimensional analysis can be used. As an exam- 
ple, let us evaluate the wavelength through the relation 
(where the 'acceleration' is proportional to the 'force' in 
the mechanical language) 



.4 
A 2 



1 



X 



(57) 



In an analogous way, we can evaluate the action J ~ 
X(A/X) 2 ~ A~ _1 , which is slighty more complicated, be- 
cause the asymptotic contribution vanishes for a > 2: in 
this case the finite, constant contribution coming from 
the motion within the 'close region' dominates. There- 
fore, J ~ A~ _1 if a < 2 and J ~ 1 if a > 2. As for the 
quantity /, appearing in the expression l|36|) for D, con- 
served case, we simply have / ~ A^4 2 ~ A = +1 . Finally, 
B(A) - Xi~ 2 . 



The nonconserved case 



From Eq. we find 

x 2 x^- 2 



D 

D - 



Xi^X 1 - 
X 2 Xi- 2 



A 1 



a < 2 
a > 2 



(58a) 
(58b) 



The phase diffusion coefficient is constant for a smaller 
than two, which means that X(t) ~ t%. For a larger than 
two, A» _1 ~ X 2 /t gives X(t) ~ t 3 <*- 2 . We can sum up 
our results, X(t) ~ t n , with 



| a < 2 
a > 2 



(59a) 
(59b) 



3a- 2 

The coarsening exponent varies with continuity from 
n = 1/2 for a < 2 to n = 1/3 for a — > oo. These results 
confirm what had already been found by one of us with 
a different approach |18| . 

2. The conserved case 
From Eq. (|36|l we have 



D 



and therefore 



A 2 A° 2 1_ 
Al+iA 1 -* ~ A 2 

A(i)~ii 



(60) 



(61) 
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The constant coarsening exponent n = 1/4 clashes 
with numerical results found in Ref. n = 1/4 for 
a < 2 and n = a/(5a — 2) for a > 2. The opinion of 
the authors of Ref. is that for a > 2 a crossover 
should exist from n = a/ (5a— 2) to n = 1/4, the correct 
asymptotic exponent. Details and supporting arguments 
will be given elsewhere [l9|. 



D. Conserved models for crystal growth 

It is interesting to consider a model of physical interest 
which belongs to the class of the full generalized Cahn- 
Hilliard equations, meaning that all the functions B(u), 
C(u) and G(u) appearing in l|29|l are not trivial. The 
starting point is Eq. iH^Hjl . 

8 t z = — d x {B(m) + G(m)d x [C(m)d x m}} 

which describes the meandering of a step, or — more 
generally — the meandering of a train of steps moving in 
phase, z is the local displacement of the step with re- 
spect to its straight position and to — d x z is the local 
slope of the step. 

We do not give details here about the origin of the 
previous equation, which is presented in but just 
write the explicit form of the functions B, G and C: 



B(m) = 



1 + m- 



-, G(m) = 



1 + (3\J\ + m 2 
(1 + (3)(l + m 2 )' 



(62a) 



ri \ l + c(l + m 2 )(l + 2m 2 ) 
C(m) = (l + c )(l + m3)3/ 2 ^ 



and define the meaning of the two adimensional, positive 
parameters appearing there: (3 is the relative strength 
between the two relaxing mechanisms, line diffusion and 
terrace diffusion; c is a measure of the elastic coupling 
between steps. 

If we pass to the new variable u = J ds C(s), we get 
Eq. 123, 



d t u = ~C(u)8 xx [B(u) + G{u)u xx ] 



(29 



whose steady states, for high-symmetry steps, are given 
by the equation u xx — —B(u)/G(u). In App.[U]we study 
the potential V(u) = J du[B(u)/G(u)] and the dynami- 
cal scenarios emerging from 8a A. We give here the results 
only. 

If c = (see Fig. 121, 8a A < 0, while there is asymptotic 
coarsening if c > (see Figs. l4l5|) . Asymptotic coarsening 
means that 8aX > for large enough A: according to 
the values of c and (3, X(A) may be always increasing 
or it may have a minimum followed by 8a A > 0: the 
distinction between the two cases is not relevant for the 
dynamics and it will not be further considered. Let us 
now determine the asymptotic behavior of all the relevant 
quantities, when c > 0. 

In the limit of large to, we have C (m) ~ m and B(m) ~ 
\jm. As for G, G(m) ~ 1/m 2 and G(m) ~ 1/m, for 



no coarsening 




FIG. 3: \{ A) for E q. lE&l . c = and = (see Eq. O). 
V(u) = — vl — u' 2 (inset) and X'(A) < (main, full circles). 
The dashed line is the harmonic potential ti 2 /2. If c = 0, 
there is no coarsening whatever is (3. 



(3 — and (3^0 respectively. Since u ~ f m ds C(s) ~ 
to 2 , m ~ \/u so that C(u) ~ y/u, B(u) ~ l/v 7 ^ an d 
G[u) - 1/m (ft = 0) or G{u) - 1/V5 (/? ^ 0). The po- 
tential V(u) has the form V(u) — u 3/2 (/3 = 0, Fig.HJ or 
V(u) ~ u (/3 ^ 0, Fig.EJ. The wavelength (see App. HJ 
is A - A 1 / 4 for /? = (Fig. gjl and A ~ y/A for /? ^ 
(Fig. |SJ). Similar and straightforward relations can be 
determined and the following general expression for the 
phase diffusion coefficient is established, 



D 



AG (A) 
A 2 



A 2 
t 



(63) 



and the coarsening exponent is finally found to be 

A(t) ~ ti (3 = (64a) 



X(t) ~ /3 7^ 



(64b) 



These results agree with both the numerical solution 
and the heuristic arguments presented in jl2j . 



E. Discussion 

In this section we have applied the results of Sec.lrflto 
find the coarsening law A(i) for some classes of models 
displaying asymptotic coarsening, that is to say having 
a negative phase diffusion coefficient D for 'large' ampli- 
tude (if the amplitude A has an upper limit A max , 'large' 
means A — > A max ). In particular, in llll AIM CI we have 
considered models entering into the classes 

d t u = B(u) + u xx =» t ~ (65a) 
5 tU = -9 ra [B( u ) + Ura ] t^^^ (65b) 

where we have also indicated on the right the relations 
leading to the coarsening laws and deriving from \D\ ~ 
X 2 /t, with D given in Eqs. I|27|) and (|36|l . 
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Sec. lIIIDl has been devoted to a class of conserved mod- 
els which are relevant for the physical problem of a grow- 
ing crystal. In that case the full expression (|35|> for D 
must be considered (the result is reported in Eqs. (|04al 
I64bll ). It is remarkable that for all the models we have 
considered, we found n < 1/4 and n < 1/2 for conserved 
and nonconserved models, respectively. It would be inter- 
esting to understand how general these inequalities are. 9 

IV. SWIFT-HOHENBERG EQUATION AND 
COARSENING 



F IG. 4: X (A) for Eq 
Vl + m 2 (2m 2 



c = oo and /3 
l)/3, with m 2 



0. 



Inset: 
l)/2. 



V(u) 
V{u) 



(VI + 4u 2 

\u\ 3 ^ 2 for large u and A'(A) > (main, full circles). For large 
amplitude, A ~ A 1 ' 4 (main, full line). The dashed line in 
the inset is the harmonic potential, u 2 /2. If c > there is 
coarsening, whose law depends on /3. For /3 = 0, X(t) ~ t 1 ^ 6 . 




FIG. 5: X(A) for Eq. c = oo an d = 00 . Inset: V(u) = 

m 2 - (1/2) ln(l + m 2 ), with m 2 = (Vl + 4u 2 - l)/2. V(u) w 
\u\ for large if and X'(A) > (main, full circles). For large 
amplitude, A w vC4 (main, full line). The dashed line in the 
inset is the harmonic potential, u 2 /2. For f3 > 0, A(f) w f 1 ^ 4 . 



Let us start from the standard Swift-Hohcnberg equa- 
tion (|57|) . 



S^u = -d*u - 2d 2 x u - (1 - 5)u - u 3 

whose linear dispersion curve is ui(q) = S — (q 2 
phase diffusion coefficient (see Eq. H44[l ) is 



(66) 
l) 2 . The 



(67) 



and it should be compared to the limiting expression, 
valid for 6 — > 0, 



D = 2q 9jMA 



(ul)) 



= - — (axS) (68) 
Zlu aq 



where the subscript E refers to the Eckhaus instability 
and prime denotes derivative with respect to q. 



Steady states satisfy the equation Csh[u] 



0, 



where Csh[u] means the linear part of the operator ap- 
pearing on the right hand side of Eq. 1)66(1 . In the one- 
harmonic approximation, u(x) = A(q) cos(qx) and the 
steady condition writes u>(q)u = u 3 implying 



iu(q) = (u 2 ) = l/2A 2 (q). 



(69) 



Passing from the standard GL/CH models fScc. lIII 
to models where V(u) have (non quadratic) maxima at 
hnite u (Sec. 



IIII B|) and to models where V{u) has no 
maxima at all at finite u (Sec. IIII C^ . the coarsening ex- 
ponents change with continuity from n = (logarithmic 
coarsening) to n — 1/2 for the nonconserved models and 
from n = to n = 1/4 for the conserved models. 

The conservation law, as expected, always slows down 
the coarsening process. Formally, this corresponds to re- 
place the action J with the quantity / in the denominator 
of D. In most cases, J is a constant while / increases as 
A K , with k > 1: a smaller D implies a lower coarsening. 
We remark that only in a very special case (models with- 
out uniform stable solutions and a < 2), I/J ~ A 2 : when 
this happens, the double derivative d xx — which charac- 
terizes the conserved models — is equivalent (as for the 
coarsening law) to the factor 1/A 2 . We stress again that 
this is an exception, it is not the rule. 



In the same spirit we get (u 2 ) = q A /2 = q~u and 
( u xx) = 9 4 ^4 2 /2 = q 4 uj, so that D reads 



D 



2 d 4 
— — [ui{q 
qu> aq 



V)]. 



(70) 



Close to the threshold, S — > 0, q = 1+K (with K < 1), 
uj(q) = 5 - AK 2 and u)' = -8K. Finally, we get 



D 



oj aq 



D f 



(71) 



so that the expression for D is equal to the well known 
expression for De- It must be noted, see Eq. Q69JI . that A 



9 The condition n < 1/2 for the nonconserved models is equivalent 
to say that |-D(A)| does not diverge with increasing A, which 
seems to be a fairly reasonable condition. 
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FIG. 6: X(A) for the Swift-Hohenberg eq. JUJl, with 5 = 1. 
The amplitude A is defined as A = [<f dx u 2 (x)/X] 1/2 . 



vanishes at q — q 2 and at q — q\ , as u> does, so that A un- 
dergoes a fold singularity at the center of the band. This 
is imposed by symmetry, since in the vicinity of threshold 
the band of active modes is symmetric. In this case we 
are in the situation with the dispersion relation • The 
phase diffusion coefficient must also be symmetric with 
respect to the center, and therefore it can change sign at 
the fold, due to this symmetry. Thus D does not have 
the sign of A' in this case. Nonetheless, in the vicinity of 
52 the sign of D is still given by A', as shown here below. 
Our speculation is that the existence of a fold is likely to 
destroy the simple link between D and A' . 

A meaningful expression for D can be found also for 

finite 5, close to 172 = \A + \/8 (let us remind that 
^(91,2) — and <72 > <Zi)- In this limit we get Lu'(q) 
u/(g 2 ) = 4q 2 



4g| and 



D = 2(4 - q 2 ] 



1 dio 
uj dq 



= D f 



(72) 



so that D is equal to a positive quantity times (Lo/dq. 
Since uj(q) = A 2 (q)/2, the sign of D is equal to the sign 
of dA/dq. This result follows from a perturbative scheme 
where use has been made of the fact that u ~ cos(qx). 
This is legitimate as long as one considers small devia- 
tions from the threshold. If 8 is not small, or if 8 = 1 we 
fall in the dispersion relation Q. As one deviates from 
c?2 towards the center of the band, higher and higher har- 
monics become active, and one should in general find nu- 
merically the steady state solutions in order to ascertain 
whether D is positive or negative. In the general case, 
we have not been able to establish a link between D and 
the slope of the steady state branch as done in the previ- 
ous sections. Our belief, on which some evidences will be 
reported on in the future, is that depending on the class 
of equations, it is not always the slope of the steady state 
solution that provides direct information on the nonlinear 
dynamics, but somewhat a bit more abstract quantities, 
as we have found, for example, by investigating the KS 
equation, another question on which we hope to report 
in the near future [2(|. Numerical solutions of the SH 
equation in the limit 8=1 reveal a fold singularity in 



the branch X(A), as shown in Fig. 



V. EQUATIONS WITH A POTENTIAL 



Some of the equations discussed in Sec. Ill Bl are deriv- 
able from a potential: it is therefore possible to define 
a function T which is minimized by the dynamics. This 
is always the case for the generalized Ginzburg-Landau 
equation l|ltj|) . which can be written as 



d t u = B(u) + G(u)u xx 



-G(u 



ST 



T[u] 



dx[\{u x ) 2 -V{u)] 



(73) 
(74) 



where V(u) is the potential entering in the study of the 
stationary solutions, i.e. V(u) = B(u)/G(u). If we eval- 
uate the time derivative of T we find 



dT 
~dt 



w 5:F 
ax-r—Ut 

8u 



dxG(u) 



5JF_ 

S u 



<0, 



(75) 



if G(u) > 0. 

The generalized Cahn-Hilliard equation l|29|l can al- 
ways be written as 



d t u = -C{u)d xx [B(u) + G(u)u xx ] = C(u)d xx [G{u) 



If C{u) = G{u), we find 
dT 



dt 



= — / dx 



d x C(u) 



8JF_ 
8 u 



< 0. 



5F 

8u 
(76) 



(77) 



We now want to evaluate T for the steady states. The 
pseudo free-energy T is nothing but the integral of the 
Lagrangian function C for the mechanical analogy defin- 
ing the stationary solutions. If E = (u x ) 2 /2 + V(u) is the 
'energy' in the mechanical analogy and J is the action, 



(78) 



where £ is the length of the one dimensional (Id) inter- 
face. We have also made explicit the dependence on A of 
the stationary solutions, u\(x). We want to determine 
the dependence on A of the pseudo free-energy, 



dT 
~d\ 



d\ 



E 



-4<o 



(79) 



where we have used that Be J = A. This result shows 
that stationary solutions with a greater wavelength al- 
ways have a lower energy. Since dynamics minimizes T, 
this result supports the idea that coarsening occurs when 
d A X > 0. 

If the curve A(^4) has extrema, there are different 
steady states with the same value of A: let us consider 
two of such states, separated by only one extremum, 



12 



X(Ai) = \{A2) = Xq. We ask what state has the lower 
free energy. If T± t 2 = J~\ u x{Ai 2 )] an d we label the differ- 
ent steady states with E = V(A) rather than with A, we 
get 



J{E 2 ) - J(Ei) 
A 

rE 2 



- (E 2 - E x ) 



1 f^ 2 

= — / dEX(E) - (E 2 - Ex) 
Ao Je 1 



E 2 — E\ 
Ao 



(A - A ) 



(80) 



where A is the average A in the interval (Ei, E 2 ). 

Therefore, if X(A) has a maximum between Ai and 
A 2 , J- 1 < T 2 , if X(A) has a minimum between Ai and 
A 2 , T 2 < T\. So, if two steady states have the same 
wavelength, the state with lower free-energy is always 
the state corresponding to a positive value of 8aX. 

Let us resume what we found in this section. For the 
gGL equations and a class of the gCH equations (C (u) = 
G(u)), there is a functional T which is minimized by the 
dynamics, dT/dt < 0; if u\(x) is a stationary solution of 
period A and we evaluate the free-energy T[u\(x)] for the 
steady branch X(A), then (d/ ' dX)T[u\(x)\ < 0, i.e. the 
free-energy decreases with increasing the wavelength. We 
are now interested to study the variation of T when the 
amplitude of u\(x) is changed, keeping the period fixed. 

Since u\(x) minimizes T, it is necessary to go to the 
second order. If u(x) — u\(x) + x(x), we get 

AT = T[u x (x) + X (x)} - T[u x (x)] (81) 

' r dx[V"(u x )x 2 (x)+ X (x) X "(x)} 

We are interested in fluctuations of the amplitude, so 
we take X ( x ) = tu\{x) and we ask if AT is an increasing 
or decreasing function of e: 



AT = - / dx[u x (x)V'{u x ) - ul{x)V"(u x )\ (82) 



What is relevant is the sign of the integral, 
Q = [ dx[u{x)V'{u) - u 2 {x)V"{u)] 



(83) 



where we have removed the subscript A to lighten nota- 
tion, keeping in mind that u{x) is a stationary solution 
of period A. A positive (negative) Q means that u(x) is 
stable (unstable) against fluctuations of the amplitude: 
we want to relate such property to the behavior of the 
curve X(A). 

Since V(u) is an even function, we can define V(u 2 ) = 
V(u). Since V = 2uV' and V" = 2V' + Au 2 V" , we find 



Q = -8 / dxu 4 V"(u 2 ). 



(84) 



In App. El Eqs. HA5IA6I) . we show that the sign of 
X'(A) is related to the sign of V"(u 2 ). For convex or 



concave potentials V(u 2 ), we can conclude that the sign 
of AT is equal to the sign of X'(A): if X'(A) > we have 
phase instability (i.e. coarsening), but amplitude stabil- 
ity; if A' (A) < we have phase stability (i.e. constant 
wavelength), but amplitude instability. If the curvature 
of V(u 2 ) changes sign, it is no more possible to establish 
such a strict relation between the signes of AT and A' (A) 
for any value of A. 



VI. SUMMARY AND PERSPECTIVES 
A. Results 

The two major results of the present work are (i) the 
derivation of a criterion for coarsening based on the be- 
havior of the steady state solutions. This criterion holds 
for several classes of nonlinear equations that are encoun- 
tered in various noncquilibrium systems. The link be- 
tween the steady state behavior and coarsening (which 
is a dynamical feature) has been made possible thanks 
to the phase diffusion equation, (ii) The exploitation of 
the phase diffusion coefficient has allowed us to derive 
the coarsening law. For all known examples which fall 
within our classes, we have captured the exact coarsen- 
ing exponent. Our analysis has allowed us to make a 
quite general statement about the law of coarsening and 
on the relevant quantities that are decisive in fixing the 
coarsening exponent. 



B. Extension to higher dimension 

Usually, an analytical derivation of the coarsening ex- 
ponent is made for some equations where their one di- 
mensional character is essential. While a link between the 
phase diffusion coefficient and the behavior of a steady 
state branch proves presently difficult to achieve beyond 
Id, the derivation of the phase diffusion equation can 
be made at arbitrary dimension. Our idea according to 
which t ~ X 2 / D is worth testing in higher dimension. If 
it works, since D contains only information on the peri- 
odic steady state solutions, it is sufficient to obtain these 
solutions to determine the coarsening law. A numer- 
ical determination of these solutions is straightforward 
and thus the behavior of D as a function of A can easily 
be extracted. Then the coarsening law can be obtained 
without resorting to a time-dependent simulation. 

It must be noted, however, that, for example, in two 
dimensions (2d), there are two independent phase diffu- 
sion coefficients. There are, besides bands, five Bravais 
lattices in 2d and one has first to determine the steady 
state solutions, which depend on the symmetries of the 
equation. If we take bands as a starting steady state 
solution, then there are two principal directions for the 
phase evolution: the phase diffusion along the band, and 
the one orthogonal to it. Usually, in the vicinity of a 
threshold (as in the case of Eq. Q) diffusion along the 
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band exhibits a zig-zag instability, while the one in the 
orthogonal direction is associated with the Eckhaus in- 
stability pj. Would only one of the two coefficients be 
relevant for the coarsening, or rather is there a compe- 
tition between the two directions? It is clear that if the 
bands coalesce by keeping their steady state like symme- 
try (bands keep their integrity), then, for the CH equa- 
tion, a logarithmic coarsening should prevail. We know, 
however, that the coarsening law in 2d is A ~ t 1/3 [2l|. 
From this, we expect that both coefficients were essential 
for the coarsening. 

One possibility would be to make use of the idea of 
diffusion law (X 2 /D ~ t) in the anisotropic case. For 
example if we have two diffusion coefficients D x and D y 
along the two directions, by absorbing these coefficients 
in x and y coordinates in the phase equation, we arrive 
to {X x /V^) 2 + (Ay/V^) 2 ~ *, where A 2 . + A 2 = A 2 , 
and D x and D y depend on A^ and Aj,. At long time 
scales, it is appealing to expect that isotropy (provided 
the starting equation enjoys the rotational symmetry) 
should be, in principle, restored and A^ ~ X y ~ t n (where 
n is the coarsening exponent), and this should complete 
the extraction of the scaling with time. It is an important 
task for future investigations to clarify this point. 

Studies on coarsening in more than Id may be fairly 
complicated, even if a lot of progress has been done in the 
comprehension of phase ordering phenomena |2l| . As a 
matter of fact, in more than Id large scale temporal sim- 
ulations which arc performed with the aim to ascertain 
the coarsening law are extremely time consuming (even 
prohibitive in several cases). Note that, as it has been 
discussed in Sec. IIII CI even in Id the cross-over to the 
true asymptotic regime may prove to be very long. Thus, 
if our idea based on the analysis of phase diffusion should 
work in 2d (and higher dimension), it would constitute 
an important way for the determination of the coarsen- 
ing exponent. We intend to investigate this matter on the 
GL equation in 2d, and then possibly on other equations. 



C. More about the phase equation 

The idea according to which the diffusion coefficient, 
even in Id, allows one to determine the coarsening law 
(although it has proven to be successful for all classes of 
equations discussed here) calls for additional comments. 
Indeed, because the phase equation dt<j> = Dd xx (f> ex- 
hibits a negative diffusion coefficient in the coarsening 
regime, the phase instability leads to an exponential in- 
crease of the phase in the course of time. One therefore 
needs to push the expansion in power series of e to higher 
orders. It is clear that the next order should lead to the 
following linear term —nd xxxx <j) with n > in order to 
prevent arbitrarily small wavelength fluctuations. Non- 
linear terms are also needed for the nonlinear saturation. 
From general considerations we expect the first nonlin- 
ear term to have the form d x (j)d xx <j>. This is dictated 
by the fact that the phase diffusion equation should en- 



joy invariance under the transformation: x — > —x and 

Why and how does the determination of the coars- 
ening exponent depend solely on the inspection of D is 
not completely understood. One possible heuristic expla- 
nation can be put forward, namely the power counting 
argument where at large scales only the first of these 
terms is decisive. It is of great importance to clarify this 
point further, for example by analyzing the renormaliza- 
tion flow of the nonlinear phase equation at large time. 

D. From coarsening to chaos 

All the evolution equations for which our general cri- 
terion based on X(A) has been derived, exhibit the fol- 
lowing dynamics: (i) they undergo perpetual coarsening, 
(ii) or they develop a pattern with a frozen wavelength 
while the amplitude increases indefinitely in the course 
of time. If, for example, the standard CH equation JfjJ is 
modified in the following manner 

d t u — -d xx [u + d xx u - u 3 ] + vud x u (85) 

where v is a parameter (taken to be positive without re- 
striction), we obtain a mixture between the CH and KS 
equations. For v — we recover the CH equation while 
for v — ► oo we obtain the KS limit (upon an appropriate 
rescaling, u — > u/v). Thus for small enough values of v a 
coarsening may be expected, while for large values chaos 
should prevail. The behavior of this solution was consid- 
ered by Golovin et al. [2^ • In the spirit of our analysis we 
can derive the phase diffusion equation, or equivalently 
we can describe the stability of steady state solutions in a 
Floquet-Bloch picture. The branch of steady state solu- 
tions in the plane (wavelength-amplitude) is monotonous 
up to v — 0.47, beyond which the branch exhibits a fold 
singularity. Numerical analysis seems to indicate that 
the transition from coarsening to non coarsening occurs 
at this value. 

While it could not be proven in general that there is 
a link between the branch of the steady state solutions 
and the instability eigenvalues, it is surprising to see that 
coarsening stops when the branch undergoes a fold singu- 
larity. On the light of the various situations encountered 
here, it is appealing to speculate that whether coarsening 
occurs or not should be related to considerations of the 
behavior of the steady state solutions only. While for the 
class of equations like l|30|) the criterion could be related 
to the slope A(A), the criterion may assume a more ab- 
stract form for more general equations. It is hoped to 
investigate this matter further in the future. 

APPENDIX A: GENERAL CONSIDERATIONS 

ON X(A) 

Stationary solutions, which play the dominant role 
in our treatment, in most cases are determined by a 
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Newton-type equation, 



we get 



F(u) 



(Al) 



where the potential V(u) = — J duF(u) is symmetric and 
has a (quadratic) minimum in u = 0. We are interested 
in periodic solutions u (x), which correspond to oscilla- 
tions of period A and amplitude A within the potential 
well V(u). We want to give general criteria to understand 
if A is an increasing or a decreasing function of A. More 
precisely, we are showing how is possible to determine 
in an easy way the behavior of A ( A) for small and large 
amplitude A. This does not cover all possibilities, be- 
cause we may find funny potentials V(u) which produce 
an oscillating function \{A): in these cases a numerical 
analysis is necessary each time. In all the other cases, 
when A (A) has no more than one extremum, the follow- 
ing analysis applies. At the end of this section we are 
also providing an exact expression for dX/dA, whatever 
is the potential V{u). 

It is trivial that in the harmonic approximation, 
V(u) ~ u 2 , A is constant: from a dynamical point of 
view, this corresponds to the linear regime. The behav- 
ior at small amplitude A is therefore denned by the sign 
of the quartic term in the potential, the third one being 



absent because V(u) is symmetric: V(u) 



■CL4U . It 



is easy to understand and straightforward to show (via 
perturbation theory j^) that 8a A has the sign opposite 
to 04: A is an increasing (decreasing) function of A if 04 
is negative (positive). 

As a general rule, if the potential is steeper than a 
parabola, A decreases with increasing the amplitude; in 
the opposite case, &aX > 0. Therefore, if V(u) w \u\ k for 
large u, A increases with A if k < 2. Using the law of 
mechanical similarity |24|. i.e. scale analysis, we can find 
that 



V(u) 



A 
A^ 



.4 



fe-i 



\~ A 



(A2) 



If V(u) increases faster than a power law, e.g. expo- 
nentially, <9aA < of course. The same is true if V(u) 
diverges at finite amplitude, e.g. V(u) = 1/(1 — u 2 ), or if 
V(u) goes to a finite value for finite amplitude, but with 
a diverging force, e.g. V(u) — —y/l — u 2 . 

If V(u) increases slower than a power law, e.g. loga- 
rithmically, or it goes to a constant for infinite u, or it 
has a maximum at finite amplitude, in all these cases A 
is an increasing function of A. 

Now let us turn to a more rigorous analysis of X'(A). 
The exact expression for the period A is 



A = 2V2 



du 



y/V(A) - V(u) 



(A3) 



If we use the formula 

_d_ " A 
dA 



duf(u, A) 



(A4) 



= j J du[Ad A f(u,A)+ud u f(u,A) + f(u,A)} 



dX V2 f A R(A)-R(u) 

dA - A J, du [V{A)-V{uW/ 2 (A5) 

where R(u) — 2V(u) — uV'(u). If, following Sec. we 
introduce the function V{u 2 ) = V(u), we find that 



R'(u) = ~4u 3 V"(u 2 ). 



(A6) 



Therefore, if V(u 2 ) is a convex (concave) function, R(u) 
decreases (increases) with u and A is a decreasing (in- 
creasing) function of the amplitude A. 



APPENDIX B: gCH EQUATION: THE SIGN OF D 

In Sec. Ill B 21 we found that for the generalized Cahn- 
Hilliard equation, the phase diffusion coefficient has the 
form D = q 2 d q ((. . . ))/-C>2- In order to establish the con- 
nection between coarsening and sign of X'(A), it is nec- 
essary to show that D2 > 0. Let us recall the following: 



D 2 = - 

d^v = 
C(u ) 



C(u )/ 
= C(-u ) > 



G(u Q ) 



(Bl) 
(B2) 

G(-u ) > (B3) 



Since C(uq) is an arbitrary even and positive function 
and v does not depend on it, it is necessary to prove 
that vd<t ,uo < 0, or, equivalently, that vd^v < 0. From 
Eq. (|B2|I . v — v (cj)) + ci + C24>, where v is a 2tt— periodic 
function: we impose C2 = in order to maintain this 
property, while c\ will be fixed later on. 

ua((j)) is the periodic, bounded trajectory of a particle 
moving in a symmetric potential well, with V(u) being 
an increasing function between u — and u = A, the 
amplitude. uq(4>) can be choosen as an odd function 
vanishing in ip — 0,±7r, with uo((f>) and —d^uo having 
the same sign as <j). If we take the derivative of uq, use the 
relation i|B2fl and integrate d^v, we recognize that d^v 
has the same properties as ua{<f). In a similar way, we 
can say that d<p<pv and d^uo are even functions vanishing 
in ±(/)o (the extrema of uq and d^v), that d^v{4>) > 
for < cj) and d^v^) < outside. 

Let us now integrate d^v: v(<p) is an even function 
having maxima in ±7r and a minimum in (f> = 0. We 
fix the constant c\ in such a way that v(<p) has zeros in 
±0o. Therefore, since v(4>) < for \(j)\ < <f) and v(<p) > 
outside, 



< 0, 



(B4) 



as we should prove. 
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APPENDIX C: THE DIFFERENT SCENARIOS 
FOR THE 'CRYSTAL-GROWTH' EQUATION 



We consider here Eq. Q30(l 

d t z = -d x {B(m) + G(m)d x [C(m)d x m]} 



(CI) 



where m = d x z and the functions B,G,C are given by 
Eq. (H3 



B(m) = 



m , 1 + /3vl + m 2 

2 ) G ( m ) = 



1 + m 2 ' ~ v '~' (1 + /3)(1 +to 2 )' 

, l + c(l + m 2 )(l + 2m 2 ) 
C(m) - 



(l + c)(l + m 2 ) 3 / 2 
If we perform a change of variable, 

m l + c(l + m 2 ) 



dsC(s) 



1 



VT 



(C2) 



we get what we named the full generalized Cahn-Hilliard 
equation 



d t u = -C{u)dxx[B(u) + G(u)u x 



(C3) 



whose steady states correspond to the trajectories of a 
particle moving in the potential 



V(u) 



du 



B(u) 

G( u y 



(C4) 



In the general case, it is not possible to invert 
the function u(m) so as to get an explicit form for 
B(u),G(u),C(u) and therefore V(u). However, it will 
be sufficient to consider the limiting expressions for V(u) 
(small and large it) in order to discriminate between the 
possible dynamical scenarios. 

Let us start by considering the case c = 0, because 
analytics is simpler. The function u(m) can be inverted, 
m = u/yl — u 2 , showing that the limit of large ampli- 
tude m corresponds for the new variable to u — > ±1. The 
potential V(u) has the form 



V{u) = (! + /?) I du 



[3 + VI -u 2 



(C5) 



where the parameter f3 > is easily recognized to be 
irrelevant, so we can assume /3 = and get the approxi- 
mate potential V(u) ~ — VI — u 2 (see Fig. [3J inset): this 
case had already been studied by the authors ^Jj , finding 
that X(A) is a decreasing function (see also the remarks 
in App.[SJ. We can conclude that no coarsening appears 
if c = 0. 

If c > 0, u(m) ~ cm 2 /(l + c) for m 3> 1 and large 
amplitude now means large u. In this limit, m ~ ^/u and 
V(u) - (1 + /3)V"/(1 + so that V(u) - y/u for 

P = and V(u) ~ (1 + /?)//? for /3 > 0. In both cases- 
potential growing as the square root of the amplitude 
or going to a constant — A ( A) is an increasing function 
and there will be asymptotic coarsening. The different 
behavior of V(u) for j3 — and (3 > suggests different 
coarsening exponents, as are actually found. 

In the previous paragraph we have considered the limit 
of large amplitude: now we are considering the oppo- 
site limit of small amplitude. When m — > it is triv- 
ial to check that u — m, but we need the third or- 
der correction, which is easy to calculate from Eq. I|C2|I : 
m ~ u + (a/2)u 3 , with a = (1 — c)/(l + c). For small m, 



B(m) 
G(m) 



(l + /3)m 



(l + /3)m 



/?Vl + m 2 
I 

2(l + /?) ? 



1 + P(l + m 2 /2) 



(C6) 



If now we pass to the variable u, 



du 



(u+V) 



2(1 + P) 



u 

T 



2(1 + /3) 



(C7) 



we find that the behavior of the curve \{A) at small A 
is fixed by the sign of the quantity a — /3/(2(l + /?)): A is 
an increasing function of A (at small A) if that quantity 
is negative, i.e. if a < /?/ (2(1 + /?)). Using the expression 
a = (1 — c)/(l + c) we finally find that c^tA > at small 
A, if c> 1/(1 + 2/3). 
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